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Looking for Signs of Anisotropic Cosmological Expansion in the 

High-z Supernova Data 
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Several problematical epochs in cosmology, including the recent period of structure formation 
(and acceleration), require us to understand cosmic evolution during times when the basis of FRW 
expansion, the cosmological principle, does not completely hold true. We consider that the break- 
down of isotropy and homogeneity at such times may be an important key towards understanding 
cosmic evolution. To study this, we examine fluctuations in the high-z supernova data to search for 
signs of large-scale anisotropy in the Hubble expansion. Using a cosmological-model-independent 
statistical analysis, we find mild evidence of real anisotropy in various circumstances. We con- 
sider the significance of these results, and the importance of further searches for violations of the 
cosmological principle. 



I. INTRODUCTION AND RESEARCH 
RATIONALE 

After a wave of successes in recent years, cosmol- 
ogy today is facing a series of potentially major turn- 
ing points. The accepted standard model of the field, 
the concordance model [3| , is a patchwork of detailed 
and observationally proven theories, interlaced with 
paradigms that are equally well accepted, though less 
well proven and far less detailed in terms of precise 
physical models. In the near future, new data will 
either bring these paradigms dramatically into focus, 
or will force great changes into the tightly interwoven 
concordance model. 

Besides the unsolved problem of baryogenesis, the 
two most problematical epochs in cosmology both 
involve similar circumstances: they seek to under- 
stand the Friedmann-Robertson- Walker (FRW) ex- 
pansion of the universe during phases in which the 
fundamental conceptual basis of FRW expansion, the 
Cosmological Principle (i.e., isotropy and homogene- 
ity) does not completely hold true. The early 
such epoch (the pre-homogenization period) is usu- 
ally handled with the Inflation paradigm [2], and 
the late such epoch (post-CMB formation of large 
scale structure) is currently addressed with paradigms 
such as Vacuum Energy 0], Quintessence 0] (e.g.. 
Tracker Quintessence (Jj), etc. Correct or not, these 
paradigms are not yet well constrained or proven in 
detail; and we take the position that it is no coin- 
cidence that the two least well understood epochs in 
cosmology happen to occur at times in which the Cos- 
mological Principle loses its validity. 

Demonstrating any direct theoretical link between 
the violation of isotropy and homogeneity, and solu- 
tions to specific cosmological problems, is beyond the 
scope of this current analysis; nor do we include here 
any investigation of the pre-homogenization period in 
the very early universe. What we focus upon here is 
the late epoch during which the CMB-era smoothness 
finally breaks down into clumpy structure. We seek 
to determine how seriously the Cosmological Princi- 



ple is broken in the recent universe, and the tool we 
use for this purpose is the collection of high-redshift 
Type la supernova data used for measuring the cosmic 
acceleration (e.g., [H), representing perhaps the best 
mapping of the expansion flow on very large scales. 
Measuring the extent of any irregularity in the smooth 
Hubble expansion would be a key step towards deter- 
mining its significance (if any) in our understanding 
of cosmological concordance. 

Theoretically, there remains room for new physics, 
since post-CMB evolution and the formation of large 
scale structure is by no means a "solved problem". 
Despite the strengths of linear gravitational collapse 
models in a Cold Dark Matter/Dark Energy (Acdm) 
universe there remain many loose ends, such as: 
the cuspy CDM halo problem, and the possible dearth 
of identifiable satellite-galaxy-sized structures in the 
local universe 01; the overabundance of high- 2 clus- 
ters predicted by low-density Acdm models, implying 
actual values of f2M higher than those in the con- 
cordance model d, 11] ; the unex pec ted existence of 
very massive galaxies at high-z [l3|; the troubling, 
persistent tendency of Type la SNe data to indicate 
best-fit values of iltot > 1 [H] and/or (for Dark En- 
ergy) w < ~1; the unknown behavior and composi- 
tion of the Dark Energy itself (not to mention the 
Dark Matter) [l2|; the possibility of non-Gaussian 
and hemisphere-asymmetric behavior in the WMAP 
data [l^ ; and the question of what is causing the lack 
of low-order multipole power in the CMB, as well as 
its rolling spectral index of fluctuations 0] . 

Searching for meaningful anisotropy in the cos- 
mic expansion is not without empirical justification. 
Well documented are the large discrepancies that have 
been historically found between different measure- 
ments of the Hubble Constant, yielding a broad and 
non-Gaussian distribution of results [14 ]: and though 
the error bars on Hq have been greatly reduced over 
the years, even results quoted as demonstrating con- 
cordance (e.g., the agreement between values for Hq 
found by WMAP and by the Hubble Key Project [1]), 
remain tempered by the fact that different types of Hq 
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measurements (e.g., combining S-Z effect and cluster 
X-ray flux measurements) stifl give somewhat discor- 
dant values [6|]. In addition, in the high-z supernova 
data used to prove cosmic acceleratioii, though the ac- 
celerating trend is statistically strong [5| , there remain 
enormous fluctuations for the individual SNe scattered 
about that trend (to be shown below). While the bulk 
of such discrepancies are doubtless due to a number 
of factors unrelated to the expansion itself - e.g., sys- 
tematic errors (theoretical and observational) , poorly- 
understood physical processes causing variations in 
SN luminosity measurement difficulties leading 
to very large statistical errors, scarcity of data, etc. 
- it remains possible that at least some of this SN 
scatter (and some of the disagreement between dif- 
ferent Hq measurements) represents real physics, and 
appears due to unknown dependencies of the expan- 
sion rate on angular position in the sky, which has 
been virtually ignored up to now by cosmological re- 
search because of the assumption of isotropy. 

Despite the broad assumption (and validity) of the 
Cosmological Principle in general, we are not the 
first to explore the possibility of meaningful cosmic 
anisotropy on the largest observable scales. Zehavi 
et al. [16| searched the early high-z supernova data 
(with marginal, positive results) for possible evidence 
that we might live in a local "bubble" of faster Hubble 
flow within a shell bounded by the local Great Walls. 
Alternatively, a detailed examination of the question 
by Lahav [T3| led to some evidence of overall homo- 
geneity, such as: agreement in shape and amplitude 
(except at the largest scales) between the 2dFGRS 
spectrum and a Acdm linear-regime perturbations 
model [3]; the near-convergence between the CMB 
dipole and the IRAS galaxy clustering dipole [13, [3 ! 
some evidence (though conflicting) of an isotropic 
distribution of very distant radio sources; an upper 
limit (produced using substantial theoretical interpre- 
tation) on anisotopy of sources contributing to the 
X-ray background; the apparent absence of big voids 
in the Lyman-a forest (covering 1.8 < z < 4, roughly); 
and, anisotropy constraints by Kolatt and Lahav 20] 
from Type la SNe. (Somewhat curious is this last re- 
sult, in which Kolatt and Lahav interpret a mild rejec- 
tion of isotropy at the 70% — 80% confldence level 
as evidence for FRW behavior at these scales.) In 
short, there are a number of lines of evidence pointing 
towards isotropy and homogeneity on very large (sub- 
CMB) scales, though none of them are completely con- 
vincing at this time. The strongest reason for believ- 
ing in extended FRW behavior long after the CMB 
epoch is still mainly a drive for concordance - a the- 
oretical motivation, not empirical proof. 

Direct probes of cosmic structure made by mapping 
the universe still remain inconclusive in demonstrat- 
ing large-scale homogeneity. Despite perennial expec- 
tations of reaching scales large enough for which struc- 
ture finally gives way to smoothness, evidence contin- 



ues to be found for apparently real structure at ever- 
increasing scales. Examples include: the large "Lo- 
cal Hole" , a significant deficit of galaxies in the APM 
survey area, re- verified with 2MASS data, and imply- 
ing possible non-Gaussian clustering on scales up to 
- 300/1-1 ^p(, p^j. _ 200 - 300 Mpc sized structure 
detected in the 2dF QSO redshift survey [l^ ; and the 
SDSS detection of the gigantic "Sloan Great Wall" , a 
structure 80% larger than the CfA Great Wall (450 
Mpc wide) in comoving coordinates jl^ . 

The implications of such results are debated: 
Mueller and Maulbetsch [23| claim that in contrast to 
other, earlier results, the supercluster and void struc- 
ture in the SDSS data is well reproduced by high- 
resolution Acdm simulations; Miller et al. (2^ demon- 
strate (using a model requiring some bias of QSO 
distribution with respect to the Dark Matter distri- 
bution) that their 2dF QSO detection of very large 
structure does not provide any evidence of collapsed, 
non-linear structures on scales larger than 100 Mpc; 
and Gott et al. [23j claim (counterintuitively, it may 
seem) that the detection of the unprecedentedly large 
Sloan Great Wall provides support for the expected 
approach to large-scale homogeneity (based upon an a 
posteriori argument that the size of Sloan Great Wall, 
despite being much larger than any previously known 
"single" structure, was not the largest possible struc- 
ture that could have fit in the SDSS survey). Never- 
theless, the overall lesson still remains: the transition 
from CMB-era smoothness to the recent, clumpy uni- 
verse is still very poorly understood, and the potential 
importance of anisotropy and inhomogeneity in the 
"Dark-Energy-dominated" epoch of cosmic evolution 
continues to be an open question. 

The immediate impetus for this analysis is the pub- 
lication of a combined, standardized list of 230 Type 
la Supernovae [2^, containing extinction, redshift, lu- 
minosity distance and angular sky position data. Sim- 
ilarly to Kolatt and Lahav [20] , we search the SN data 
for statistical evidence of a lack of uniformity in the 
Hubble flow. What is different, however, is that we do 
not interpret our results according to any particular 
cosmological model. While they express their results 
in terms of variations in Hq, Qm, etc., we instead 
ask a bare statistical question: "Have some regions 
of the universe been expanding faster than others?" , 
once the average Hubble evolution with respect to z 
is (empirically) removed. This has two advantages: it 
keeps our analysis much more independent of theoret- 
ical assumptions; and we avoid the problem of having 
to dilute data from a limited, noisy sample by divid- 
ing the statistical power of the results among 4 four 
different model parameters. 

One other aspect of our analysis is that we sepa- 
rately consider SNe at "intermediate-z" (.01 < z <'--^ 
.1 - .2), vs. "high-z" (z .1 - .2). We use different 
types of analyses because the former is distributed 
more evenly, giving better sky coverage; though the 
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latter is of more interest here, since the higher-z 
SNe occupy more cosmologicahy-significant volumes 
of space, and correspond to look-back times closer to 
the onset of Hubble-flow acceleration. 

Finally, we note that this conference proceedings 
paper is just a general overview of our results; we plan 
to present a more detailed discussion of our statistical 
methods and results in a future journal article. 

II. PREPARATION OF DATA FOR 
ANALYSIS 

The data analyzed here is taken from the work 
of Tonry et al. [25|, a compilation of Type la SN 
data from many sources, including the High-z Super- 
nova Search Team (HZT), the Supernova Cosmology 
Project (SCP), and other researchers throughout the 
years. The total sample includes 230 SNe, presented 
in as uniform a fashion as possible, with a common 
calibration. 

For most of our cosmological analysis, we apply the 
same data cuts as in [2^: we remove all SNe with 
z < 0.01 (to minimize scatter from non-Hubble-flow 
peculiar motions), and we remove SNe known to be 

I 



In most analyses, these A^'s would then be used to 
find a best-fit cosmological model to be interpreted in 
terms of 51^, ^a, etc. But what we are interested in 
here is not finding best-fit model parameters, but in 
studying the scatter of the data around the average 
cosmological expansion (whatever that is, assuming 
it to be well defined), for potential signs of cosmic 
anisotropy and inhomogeneity. To do this in a purely 
statistical, model-independent way, we perform a very 
simple fit to the data, and then subtract this fit, f{z), 
from the A^'s to compute the "modified residuals": 

S{m - A/), = 5, - [A(m - M) - /(z)], . (4) 

These J^'s are reasonably independent of z, and can 
thus be used for SN statistical analysis. 

A polynomial fit to the Si for the 172 SNe is 
shown in Figure [l] with a best-fit function of /(z) = 
(0.586z — 0.547z^). The essential physics is fairly well 
modeled (for this hmited range of z) with just a 2"^- 
order polynomial, showing both the current cosmic ac- 
celeration and the earlier deceleration. This fit yields 
a of 230; a little bit high, but good enough for our 
purposes. The cosmological fits in [25| do somewhat 
better, but they achieve that partially by artificially 



heavily extinguished (i.e.. Ay < 0.5 mag, for SNe 
with listed Ay values). These cuts reduce the sample 
set from 230 to 172, but successfully eliminate many 
dramatic (and likely spurious) outliers from the cos- 
mic expansion rate fits. Later, we will divide the data 
into groups of "mid-z" and "high-z" SNe; but for now 
we consider the entire sample. 

The distance modulus (in Mpc) for each supernova 
is given as [2g]: 

(m-Af) = 25+51og(di) = 25+5 log(dii/o)-5 log(i7o) , 

(1) 

where dL is the supernova's luminosity distance, and 
lists of \og{dLHo) for all SNe are given in [2B |. 

The theoretical luminosity distance in an empty or 
coasting (i.e., qo = 0) universe would be [26j : 



with an expression similar to Eq. [l]for calculating the 
theoretical distance modulus, (m — M)coasting- Any 
cosmic acceleration or deceleration is then found by 
computing the residuals, as follows (for the i*'' SN): 



(3) 



■1.5^ 
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FIG. 1: Residuals of A(m — M) for Supernovae, plot vs. 
best- fit 2"'*-order polynomial, y = 0.586z - 0.547^^. 



enlarging the error bars (which are somewhat heuris- 
tically generated in the first place) to include a ve- 
locity uncertainty (or dispersion) of 500 km/sec. We 
do not do this here, since such a procedure actually 
obscures the effect which we are trying to study (i.e., 
their "noise" is our data). 

For any fit with a reasonable number of parameters, 
however, the SNe will always have tremendous scatter 
about the best-fit trendline. In this case, the index of 



A(m - M), = Ai = [(m - AT) - (m - M)coasUng]t = 5 [log(dLffo) - log(cz) - log(l + -)], . 
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FIG. 2; Modified residuals, 5{m - M), from Fig.[T]fit. 

fit [13 i? is only 0.32, indicating that ~70% of the 
variation in the data, roughly speaking, has nothing 
to do with this (or any similar) fit. Tonry et al. [25| 
note the extreme noisiness of this data, and deal with 
it by binning the SNe, and taking medians over each 
redshift bin. While this may be useful for the analysis 
of cosmological models, it once again is a procedure 
which masks the physics that we wish to study here. 

A plot of our fit-removed 5iS is shown in Fig. [21 de- 
picting a similar-looking (though now visibly random) 
scatter. This is the main data that we analyze in this 
paper. 
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FIG. 3: Hammer- AitofT plot of 5{m - M) for SNe at .01 < 
z < .2 (with bubble size oc z), color-coded depending upon 
5 (blue for low values, yellow for medium values, red/pink 
for high values). 
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III. PRESENTATION OF RESULTS 

Here we give a brief overview of some results of our 
Supernova analysis. 

First, considering the "mid-z" (i.e., .01 < z <~ 
.1 — .2) SNe data, we divide the SNe into 3 groups, 
depending upon their 6i values: "high", "mid", or 
"low" . A plot of these data on the sky is shown in 
Fig. [31 and a modal analysis (using Spherical Har- 
monic modes) is depicted in Fig. [H shown along with 
the modal decompositions of 300 randomly simulated 
skies, for comparison. 

Some significant results here: the possible detec- 
tion of a moderate dipole (larger than 80% of simu- 
lated skies); no evidence of a significant quadrupole 
mode (smaller than ~ 70% of simulated skies); and 
the possible detection of one or two especially large 
anisotropy modes, particularly the Y33 mode (cosine 
phase) . The structure of this specific mode is demon- 
strated in Fig. [51 and future analysis with more data 
will be crucial in determining whether this is a real 
sign of anisotropy in the cosmic expansion, or merely 
a random event from the decomposition of this data 
set into a large number of modes. 

In an attempt to interpret the dipole mode present 
in the data (if it is real). Figure [5] compares the di- 
rection of this dipole with a variety of other known, 
cosmologically-significant dipoles or significant di- 
rections - such as the CMB Dipole, and dipole 
hotspots/coldspots from 2Mass, IRAS, Local Group 
velocity, etc.; and the Milky Way center, the Great 



FIG. 4: Spherical harmonic modal decomposition of S{m — 
M) for .01 < z < .2 SNe (pink), vs. modal decompositions 
for 300 simulated random skies (blue). 

Wall, the Great Attractor, the Perseus-Pisces Su- 
percluster, the Supergalactic Plane, etc. Simple in- 
spection does not indicate any clear, obvious align- 
ment of the dipole we have derived from the SN data 
with any of these other dipoles, though it is possi- 
ble that the Supergalactic Plane does run through 
the hot/coldspots of this apparent SN dipole. In any 
case, future supernova data should more precisely de- 
termine the direction and magnitude of the SN dipole 
(if it exists), and determine whether or not there is 
any real correlation with other cosmic dipoles. 

Now we consider the "high-z" (i.e., z > .1) SNe 
data, again dividing it up into three groups, with 
"high" , "mid" , and "low" S.i values. A plot of these 




FIG. 5: Single largest mode from the decomposition in 
Fig. m YiiCosine (real mode, cosine phase). 



1301 



FIG. 6: Comparison of dipole direction from the decom- 
position in Figs. [3]&|4l vs. various galactic, intergalactic, 



FIG. 8: Same data as in Fig. [71 re-partitioned into more 
smaller, more numerous groupings. 



IV. CONCLUSIONS 
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FIG. 7: Hammer-Aitoff plot of 5{m - M) for SNe at z > 
.1 (with bubble size cx z), color-coded depending upon 5 
(blue for low values, yellow for medium values, red/pink 
for high values). Data partitioned into 4 regions reflecting 
groupings in space and z. 



data on the sky is shown in Fig. [71 with the sky 
roughly divided into 4 quadrants, for statistical com- 
parison. (The kind of modal analysis done earlier for 
the intermediate-z SN data is not as appropriate here, 
due to the lack of any real "all-sky coverage" by the 
data.) 

Statistical tests (t-tests ^Jl) between the upper- 
left and lower-right quadrants (the quadrants with the 
most and highest- 2; SN data) show mild, positive ev- 
idence of an asymmetry between the quadrants, with 
an effective statistical significance of ^ 1 — 1.8cr (de- 
pending upon the specific data cuts used). 

Lastly, as is shown in Fig. [51 we re-partition the SNe 
data into smaller, more specific groupings (partition- 
ing done somewhat arbitrarily, but with the partitions 
being made as similar as is feasible in sky area). In 
this case, statistical tests (AN OVA tests [13]) again 
indicate mildly positive results (~ 1 — 2a, depend- 
ing upon the specific data cuts) for real differences 
between the 5i values - and thus the expansion rate 
histories - for SNe in different parts of the sky. 



We summarize these overall results as follows: 

1. Given that Friedmann Robertson- Walker Cos- 
mological Expansion depends upon the Cosmological 
Principle (Isotropy and Homogeneity), and that these 
assumptions break down in the recent universe (the 
Structure Forming and Accelerating Epoch), it is im- 
portant to test the extent of this breakdown. 

2. High-z Supernovae are likely the best probes 
for these tests, as long as angular information is con- 
sidered and analyzed, not just sky-position-averaged 
behavior as a function of z. 

3. For SNe at .01 < z < .2, we find some evidence 
of a dipole (which the Supergalactic Plane may be 
passing through), as well as a largest (real Spherical 
Harmonic) Anisotropy mode of Yaa Cosine. But both 
findings are difficult to quantify in terms of statistical 
significance. 

4. For SNe at z >^ .2, we find some mild pos- 
itive evidence for a "Dipole-like" Expansion Rate 
Anisotropy in opposite regions of the sky; and similar 
evidence for anistropies between smaller subdivided 
regions of the sky. But large gaps in sky coverage 
make these results hard to evaluate conclusively. 

5. As more and better Supernova data are obtained 
(especially more all-sky coverage), we will be able to 
place more significant statistical limits on these po- 
tential anisotropics in the cosmological expansion. 
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